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Abstract. The properties of the high-temperature phase of Yang-Mills theory in Landau gauge are inves- 
tigated by extending an earlier study on the infinite-temperature limit to finite temperatures. To this end 
the Dyson-Schwinger equations for the propagators of the gluon and the Faddeev-Popov ghost are solved 
analytically in the infrared and numerically at non-vanishing momenta. Gluons, polarized transversely with 
respect to the heat bath are found to comply with the Gribov-Zwanziger and the Kugo-Ojima scenario, 
while longitudinally polarized gluons are screened. Therefore the high-temperature phase is strongly in- 
teracting. It is furthermore conjectured that Yang-Mills theory undergoes a first-order phase transition. 
Indications are found that at high temperatures the thermodynamic properties are nearly those of an ideal 
gas, although long-range interactions prevail. 
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1 Introduction 

Many properties of the high-temperature phase of QCD 
and even of pure Yang-Mills theory are not yet under- 
stood, especially the infrared behavior and thus the fate of 
confinement. The equation of state, as found in lattice cal- 
culations, exhibits an almost trivial Stephan-Boltzmann 
behavior of the ideal gas pp. This is also found in so- 
called weak-coupling expansions 0- On the other hand, 
it has already long ago been argued that the microscopic 
processes cannot be purely perturbative [3]. This is ex- 
emplified by the non- vanishing spatial string-tension rj]. 
In addition, in the limit of infinite temperature, at least 
part of the gluon spectrum shows confinement |S] in ac- 
cordance with the scenarios of Kugo and Ojima and 
of Gribov and Zwanziger • Also at finite temperature 
indications of such effects have been seen using various 
techniques This clearly demonstrates that non- 

perturbative quantum effects persist in that limit. 

The objective here is to understand the properties of 
the gluons and Faddeev-Popov ghosts, in particular in 
the infrared, in the high-temperature phase. This work 
is therefore an extension of studies at infinite tempera- 
ture [3] and complements studies at temperatures below 
the phase transition |12j . 

To accomplish this task, the Dyson-Schwinger equa- 
tions (DSEs) 13 for the propagators of the gluon and 
the Faddeev-Popov ghost are investigated. A clear sepa- 



ration of soft and hard degrees of freedom is found. In- 
dications will be discussed that this separation of scales 
is the reason why a microscopically highly non-trivial and 
especially nonperturbative theory can exhibit quite simple 
macroscopic properties as is the case for the equation of 
state. 

The paper is organized as follows. In section [21 Yang- 
Mills theory and signals of confinement are briefly re- 
viewed. The Dyson-Schwinger equations are introduced in 
section[3] An analytical treatment of the infrared sector is 
given in section^ The consequences of the necessary trun- 
cations and renormalization will be discussed in sectional 
Numerical solutions beyond the infrared will be presented 
in section El This includes comments on the thermody- 
namic potential and Schwinger functions. A discussion of 
the results and the phase structure are finally given in 
section A summary and concluding remarks close the 
paper in section |S] Some technical issues are deferred to 
two appendices. 



2 Aspects of Yang-Mills theory 

The following investigations are restricted to pure Yang- 
Mills theory as substantial evidence exists that the non- 
perturbative features of QCD are generated in the gauge 
sector. Hence, the theory studied here is an equilibrium 
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Yang-Mills theory governed by the Euclidean Lagrangian 1 

C — - F a F a 4- r a r) D ab r b 

F^,v — d^A^, — d v A® — gf abc ' A h ^A c vl 
Df = S ab d^ + gf abc Al . 

Here F" denotes the field strength tensor, D ab the covari- 
ant derivative, g the gauge coupling and f ahc the structure 
constant of the gauge group, which is SU(3) unless oth- 
erwise noted. A" denotes the gluon field, and c a and c a 
are the Faddeev-Popov ghost fields describing part of the 
quantum fluctuations of the gluon field. We choose the 
Landau gauge, which for technical reasons is best suited 
for the purpose at hand |14j . 

The main focus of the present work is the fate of con- 
finement at temperatures above the phase transition. It 
is therefore necessary to be able to extract information 
about confinement. Part of this information is encoded 
in the pertinent 2-point functions. The corresponding sig- 
nals will be only listed here. For a brief introduction to 
confinement in covariant gauges see ref. |lfi| . 

A sufficient criterion for the presence of confinement is 
based on the fact that no Kallen-Lehmann representation 
of a particle exists if its spectral function is not positive 
semi-definite. It is then not part of the physical spectrum 
and thus confined, c.f. |17] . This is the case if the corre- 
sponding propagator D vanishes at zero momentum, 

lim D(p 2 ) = 0. (1) 

This is essentially the behavior expected for the propaga- 
tor of a confined gluon. 

A second criterion stems from possible confinement 
mechanisms. The Kugo-Ojima scenario puts forward 
the idea that all colored objects are BRST charged and 
thus unphysical. One precondition for this mechanism is 
an unbroken global color charge. In the Landau gauge this 
condition can be cast into JS| 

lim p 2 D G (p 2 ) -» oo, (2) 

where Dq is the propagator of the Faddeev-Popov ghost. 
This scenario necessarily also implies the condition (JIJ for 
the gluon propagator. 

In the Gribov-Zwanziger scenario f3, entropy argu- 
ments are employed to show the dominance of field con- 
figurations close or on the Gribov horizon in field config- 
uration space. This scenario predicts again condition i|2"jl. 
For an infrared constant ghost-gluon vertex, which is sup- 
ported by lattice calculations ^Hj and semi-perturbative 
calculations [23j , condition (JJJ follows as well for the gluon 
propagator [Sj. 

Also intuitively it is clear that a strongly divergent 
ghost propagator at zero momentum can mediate confine- 
ment. Such an infrared divergence relates to long-ranged 

1 The hermiticity assignment for the ghosts is different from 
the conventional, although valid for Landau gauge 14 . 



spatial correlations. These are stronger than the ones in- 
duced by a Coulomb force since the divergence in momen- 
tum space is stronger than that of a massless particle. 
An infrared vanishing gluon propagator is also intuitively 
linked to confinement, as p 2 — puts a non-interacting 
gluon on-shell. Thus the gluon does not propagate and is 
confined. 



3 Dyson-Schwinger equations 

The DSEs |13U14| form an infinite tower of coupled non- 
linear integral equations for the Green's functions of a 
given theory. Therefore, in general only a truncated set 
can be solved in practical calculations. In the following 
we aim at a closed set of equations for the pertinent 2- 
point functions. In Landau gauge and at finite tempera- 
ture these are the ghost propagator 

D G {p) = (3) 
and the gluon propagator |15| 

dM = PtM ^^ +PlM i M^1. (4 ) 

where Pt and Pl are projectors transverse and longitu- 
dinal w.r.t. the heat-bath. Eqs. J3J and |@J define the di- 
mensionless dressing functions G(p^,p 2 ), Z(p^,p 2 ), and 

The derivation of the Dyson-Schwinger equations is a 
straightforward, but tedious task. As in the previous in- 
vestigations |51ll2| . we follow here ref. [2] and keep only 
the equations for the propagators. Furthermore we neglect 
one-particle-irreducible two-loop diagrams and assume a 
perturbative color structure. Furthermore, the ghost-gluon 
vertex is taken to be bare, in accordance with recent in- 
vestigations ■ 1 0(1201 . The construction of the various three- 
gluon vertices will be discussed in sectional 

The equations arc obtained from the vacuum equations 
(see e.g. ref. ^3]) by application of the Matsubara formal- 
ism ^ni- To obtain scalar equations for the (infinite) set 
of Matsubara frequencies of the dressing functions Z and 
H of the gluon propagator Q , the gluon equation is con- 
tracted with the generalized projectors P?^ an d P^^, 
respectively, defined by 2 

P K„ = C p r^ + (C - i) (<V - V)M , (5) 

^ = ^ + «-l)(^ + ^). (6) 

The choice of the projectors was made so as to obtain 
a well-defined 3-dimensional limit. The parameters £ and 

2 These are chosen differently from the ones in ref. [5] to 
avoid an inconvenient kinematical singularity. In the limit of 
infinite temperature, the resulting equations for the soft modes 
are the same. 
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£ allow to vary the projection continuously in order to 
investigate the amount of gauge symmetry violations. The 
dependence on £ vanishes in the inhnite-temperature limit 



Employing the projectors J5J and I© yields the finite 
temperature DSEs as 



G(p) 



(2 



]T / dtfd |q| (A r (p, g)G(g)Z(p - g) 

/ *-l _. , * 



71— — OO 



-A L (p,q)G(q)H(p-q) 



(7) 



HG i rpHH 



H(p) 



£Z 3L +T n ^ +T 



g 2 TC A 



c>o „ 

£ I d0d\q\(P(p,q)G{q)G(p + q) 



n— — oo 



z(p) 



+N L (p, q)Z(q)Z(p + q)+ Ni(p, q)H(q)Z(p + q) 
+N 2 (p, q)H(p + q)Z(q) + N T (p, q)H(q)H(p + q) 

(8) 

= Z 3T + T GH + T GG 



oo „ 

]T / d0d\q\(R(p,q)G(q)G(p + q) 



9 2 TC A 
' ^ .. 

-M L (p, q)H(q)H(p + q)+ Mi(p, q)H(q)Z(p + q) 
-M 2 (p, q)H(p + q)Z{q) + M T (p, q)Z(q)Z(p + q) 

pS(C - 1) l ZaL i 



2p 2 



H{p) 



(9) 



Here 5 a d C A = f dbc fabc = S a d N c = <5^3 is the adjoint Casimir 
of the gauge group. The summation runs over all Matsub- 
ara frequencies go = 2-7rTn. This (truncated) set of DSEs 
is graphically displayed in Fig. ^ 

The C an d £ dependence of the DSEs is acquired by 
using J5J and JBJ). At Q = £ = 1 the original form of the 
equations is recovered. For po = 0, equation JHJ is only 
superficially dependent on £. As all integral kernels are 
in this case proportional to £, the dependence can be di- 
vided out for £ ^ 0. Then, only the implicit dependence 
through interactions with the hard modes remains. The 
latter effect vanishes in the infinite-temperature limit and 
JHJ becomes the equation for the Higgs field of 5 , inde- 
pendent of £. 

The finite-temperature theory is renormalizable. Thus, 
explicit wave-function renormalization constants Z3 , Z?,l , 
and Zyp have been introduced and will be discussed in 
section Concerning the vertex renormalization, Z\ = 1 
has been employed for the ghost-gluon vertex |22J - The 
thrce-gluon vertex renormalization constant Z\ cannot be 
calculated in the present truncation, but it is finite and set 
to 1 [251121] . The kernels A T) A L , R, M T , M u M 2 , M L , 
P, Nt, Ni, N2, and Nl contain trivial factor such as the 
integral measure. The integral kernels are quite lengthy, 



and will not be displayed here. They can be found in ref. 

The tadpoles T y are used to cancel spurious diver- 
gences in much the same way as in ref. see sectional 
However, also in the soft equations, they can possess fi- 
nite parts at finite temperature. In case of the longitudinal 
equation, these can be absorbed into the mass renormal- 
ization discussed below. In the transverse equation, these 
are at £ — 3 completely contained in the loop terms, and 
their continuation away from £ = 3 is thus arbitrary. Also 
these contributions scale at best only as 1/p 2 in the in- 
frared and thus turn out to be irrelevant. Therefore, they 
are dropped, especially as no simple prescription as the 
one given later in Eq. I|23ll can remove the related spuri- 
ous divergences 

Close inspection of the equations reveals further that 
the dressing functions can only depend on \po\ and 
The corresponding symmetry, under pq — * — po, is used to 
reduce the number of equations significantly. 

Eqs. I|7I9(I depend on the coupling constant g as the 
only parameter. In turn, g depends on the renormaliza- 
tion scale /i, which can be chosen arbitrarily at any fixed 
temperature T. It directly enters into the definition of the 
infinite-temperature limit, as the effective 3-dimensional 
coupling constant g 3 depends on g 26 . In the simplest 
case, g 2 ~ g 2 (/i)T. Comparing Eqs. (|7I9|I to those of the 
infinite-temperature limit [H], the constant of proportion- 
ality has to be 1 in this truncation scheme in order to ob- 
tain a smooth infinite-temperature limit. Hence, it only re- 
mains to choose the temperature dependence of \i. As dis- 
cussed in appendix lAl a t'Hooft-like scaling 27 with g 2 T 
fixed is employed to obtain a smooth infinite-temperature 
limit. 



4 Infrared properties 

Asymptotic freedom turns out to be advantageous for the 
infrared as well. Just at the phase transition, the n = 1 
Matsubara frequency has already an effective 'mass' po = 
2irT c w 1.7 GeV. By the Appelquist-Carazzone theorem 
[25] the hard modes are suppressed by powers of |p|/po m 
the infrared. Thus, hard modes do not contribute signifi- 
cantly. Therefore, the infrared behavior of the soft modes 
is the same as in the infinite-temperature limit and given 
by power-laws, 

G(p) = A gP - 2K , 
Z(p) = A zP - 2t , 
H{p) = A h p- 21 . 

The longitudinal function H behaves mass-like with 1 = 1, 
The exponents of the ghost and transverse gluon propa- 
gators are related as |51I29] 



1 ( 1 
— \ t+ - 

2 V 2 



(10) 



As in the infinite-temperature limit, always two solutions 
are found. One has k = 1/2 while the other has a weakly 
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Fig. 1. The truncated Dyson-Schwinger equations at finite temperature. The dotted lines denote ghosts, the dashed lines 
longitudinal gluons and the wiggly lines transverse gluons. Lines with a full dot represent self-consistent propagators and small 
dots indicate bare vertices. The open circled vertices are full and must be constructed in a given truncation scheme. A bare 
ghost-gluon vertex and slightly modified bare gluon vertices have been used here. Note that for soft modes the ghost-longitudinal 
and the 3-point coupling of three longitudinal and of one longitudinal and two transverse gluons vanish. 



C-dependent value, being k = 0.39760 at C = 1 (5123] . The 
later is the more likely solution, at least in the infinite- 
temperature limit 5 . Hence the infrared exponents are 
independent of temperature and therefore at least glu- 
ons transverse w.r.t. the heat-bath are confined above the 
phase transition. Studying each equation in detail, it turns 
out that these quite general statements are implemented 
very differently in each equation. 

In case of the hard modes, there is no purely soft con- 
tribution due to momentum conservation at the vertices. 
Thus they decouple and a self-consistent solution is that 
all hard-mode dressing functions are constant in the in- 
frared. 

In the case of the soft ghost equation Q, all contri- 
butions become constant in the infrared. Thus they can 
be canceled by the (renormalized) tree-level value. Only 
the subleading behavior remains, as in the case of the 3- 
dimcnsional and 4-dimensional theory |51l21| . As the sub- 
leading contribution of the purely soft term dominates the 
hard terms, the same behavior as in the infinite-temperature 
limit emerges. 

In case of the transverse equation the hard mode 
contributions give rise to mass-like 1/p 2 terms at best. 
These terms are subleading compared to the soft ghost- 
loop. Thus the infrared is dominated by the latter and the 
same infrared solution as in the infinite-temperature limit 
is found. 

Finally, the longitudinal equation JHJ) is quite different 
from its equivalent in the 3-dimensional or 4-dimensional 
case. In the prior case, it was dominated by its tree-level 
mass, while in the latter it is identical to the transverse 
one. In the present case, the absence of pure soft interac- 
tions with ghosts, which could generate divergent contri- 
butions as in the transverse equation, leads to dominance 
of the hard modes 3 . These provide mass-like contributions 



3 The pure soft interactions represented by the kernels Ni 
and N2 only generate a constant term in the infrared. With- 



in the infrared. Inspecting e.g. P at £ = 1 

Iql of? sin , . 

P(0,q ,q,p)= q ° — , 11 

\p\ q 2 (q 2 o + {p + q) 2 ) 

and using the fact that the hard dressing functions are 
constant in the infrared, directly leads to a mass-like be- 
havior due to the explicit l/|p| 2 factor and the finiteness 
of the remaining expression. Thus, the 3-dimensional mass 
of the Higgs sector is generated dynamically by the inter- 
action with the hard modes. 

In conclusion, the infrared exponents of all soft mode 
dressing functions do not depend on the temperature. Only 
the infrared coefficients change. 

These observations imply that, by using a finite num- 
ber of Matsubara frequencies, it will not be possible to 
obtain the infrared 4-dimensional vacuum solutions. This 
would only be possible, if the infrared limit of the high- 
temperature and the low-temperature phase would have 
been both encoded in the soft interactions alone or, as the 
electric screening mass, would already be induced by a sin- 
gle hard mode. However, any finite number of modes will 
not be able to generate a divergence stronger than mass- 
like in the longitudinal equation, equivalent to the soft 
ghost contribution in the transverse equation. Hence, the 
4-dimensional behavior can be established only by infinite 
summation. In this case, the generated mass has to diverge 
even after renormalization at p = to obtain the vac- 
uum solution, providing over-screening instead of screen- 
ing. The possibility that the ghost-gluon vertex changes 
due to temperature resulting in a coupling of soft ghost 
modes to soft longitudinal gluons at low temperature but 
not at high temperature is unlikely, since a bare ghost- 
gluon vertex is already sufficient at low temperatures and 
in the vacuum jl2U21| . 



out a tree-level mass, the soft tadpoles generate a contribution 
which can be removed by renormalization. 
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5 Truncation 

Concerning the truncation chosen, all arguments proven to 
be successful in the vacuum and the infinite-temperature 
limit [51l21U25| also apply here. This includes the require- 
ment that, due to the Gribov condition, the dressing func- 
tions G, Z, and H have to be positive semi-definite. Hence 
the discussion here will only cover the additional aspects 
due to finite temperature. The main additional trunca- 
tion is to include only the finite set [— N + 1,N — 1] of 
Matsubara frequencies. 

For large but still finite T the 3-dimensional limit is to 
be made explicit. Thus, all the problems encountered in 
the 3-dimensional theory [2] persist, including the neces- 
sity of a modified soft 3-gluon vertex. Hence, all modifi- 
cations 4 of the pure soft terms will be left as in the case 
of [5], except for the tadpoles in the Higgs equation ©. 
Due to the dynamically generated mass, it will be neces- 
sary to alter this behavior. As no further information on 
the interaction vertices involving hard modes is available, 
the corresponding vertex functions will be assumed bare. 
Due to the large mass of the hard modes, such a tree-level 
ansatz should be justified. 

The hard mode contributions induce additional spuri- 
ous divergences. As they contribute to the two-point func- 
tions mostly at mid-momenta, when viewed from the 4- 
dimcnsional perspective using p 2 = Po+p 2 , it is expected 
that in the chosen truncation, they are much harder to 
compensate. This is indeed the case. The subtraction pre- 
scriptions are listed in appendix El The technically quite 
involved construction is described in detail in ref. |25| . 

After removal of the spurious quadratic divergences, 
the resulting equations are still logarithmically divergent 
for N — > oo. The integrals themselves are convergent, but 
the Matsubara sum is not, as the terms scale like 1/qo- 
Hence the sum diverges for an infinite number of Matsub- 
ara frequencies. This is the 4d logarithmic divergence and 
thus the usual one of Yang-Mills theory , which has to 
be renormalized appropriately. To this end, three regions 
of momentum are to be distinguished: 

At p < 2ttT, all terms in the Matsubara sums behave 
essentially the same due to the Appelquist-Carrazone de- 
coupling theorem. It is valid as the hard modes behave 
essentially tree-level-like. As the contributions of the hard 
modes behave as 1/qo, the final result will depend on the 
number of Matsubara frequencies included. Thus, it is nec- 
essary to renormalize in order to be independent of the 
cutoff, which is in this case imposed by the number of 
Matsubara frequencies included. The renormalization pro- 
cedure will be implemented below. By this approach, the 
results can be made quite reliable in this regime. 



4 Recent investigations indicate that the three-gluon vertex 
may have an infrared divergence |30] . As pointed out also in 
ref. |30| . this does not affect the infrared regime, as this range 
is dominated by the ghost loop. This is likely also the case at 
finite temperature in the transverse equation. Due to the chro- 
moelectric screening it is conceivable that this behavior does 
not persist in the longitudinal equation at high temperatures. 



At 27rT < p < 2ttT(N — 1), more and more Matsubara 
terms depart from their 1/qo behavior to a 1/p behavior. 
As the external momentum p becomes large compared to 
the effective mass qo of a hard mode, the mode becomes 
dynamical and behaves like a massive 3-dimensional par- 
ticle. The results are still quite reliable in this region, since 
also for N — > oo, at any finite momentum p only a finite 
number of Matsubara frequencies are dynamical. It be- 
comes less and less reliable when approaching the upper 
limit 2ttT(N- 1). 

At p > 2ttT(N — 1), the situation changes drastically. 
Opposite to the case N = oo, all Matsubara modes are 
dynamical and their contributions will scale as I /p. Thus 
the number of Matsubara modes will now enter linearly 
instead of logarithmically. This is an artifact of cutting 
off the Matsubara sum. In this region the results are not 
reliable. However, as the sum is still finite and suppressed 
by 1/p, the contribution is subleading with respect to the 
tree-level term and thus the system of equations can still 
be closed consistently 5 . 

By renormalization, these artifacts can be reduced, if 
not completely removed at sufficiently small momenta. In 
that sense the finite Matsubara sum approximation is a 
small- (3-)momentum approximation. However, there are 
a few subtleties involved concerning the renormalization 
of a truncated Matsubara sum. 

In the vacuum case |14H23) . the DSEs can be renormal- 
ized using a momentum subtraction scheme (MOM). This 
approach fails if the large momentum asymptotic value of 
the dressing functions is a constant different from 0. In 
the case of a truncated Matsubara sum, the ultraviolet 
behavior is that of a massive 3-dimensional theory. Thus 
the self-energy contributions vanish as 1/p in the ultravi- 
olet Up- Therefore the dressing functions F = G, Z, H are 
dominated by the tree-level term 

km F(po, |p|) - (12) 

|p| — >oo Z/Q 

Here Z3 is generically the wave function renormalization 
of F. For a finite number of Matsubara frequencies Z% is fi- 
nite. Therefore the renormalization will here be performed 
by explicit counter-terms. This is discussed below. 

A second point is the mass and mass renormaliza- 
tion necessary for the soft longitudinal mode. The soft 
mode with frequency po = of the Aq component of the 
gauge field transforms homogeneously instead of inhomo- 
geneously under gauge transformations. Therefore gauge 
symmetry permits to add a term 

Sm 2 A 2 (0,p) (13) 

to the Lagrangian. In the vacuum such a term is forbid- 
den by manifest Lorentz invariance. At finite temperature 
in the Matsubara formalism, this is no longer the case, 
and such a term could in principle be present. This term 

5 A further argument is that otherwise the corresponding 
finite 3-dimensional theory would be ill-defined which is not 
the case. 
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replaces the PqAq term of the hard modes, which stems 
from the Aod^Ag term in the Lagrangian for po ^ 0. 

Concerning the counter-terms, the wave-function renor- 
malization is performed by adding the counter-term 



5Z s (A%d»d v J%-Al&A%) 



(14) 



to the Lagrangian. In the case po = 0, the first term is not 
present for the longitudinal mode Aq. Its place is taken 
by IjlHjl . This implies a relation between the wave- function 
counter-term 5Z 3 and the mass counter-term 5m 2 , which 
cannot be exploited here due to the truncation of the Mat- 
subara sum. Therefore an independent renormalization of 
the wavefunction and the mass of the soft longitudinal 
mode is necessary and will be performed. 

A last point concerns the implications for the counter- 
terms due to the truncation of the Matsubara sum. As the 
divergence structure must be the same as in the vacuum, 
the counter-terms must be the same for all frequencies. 
This is no longer the case when the Matsubara sum is 
truncated, as can be seen directly by counting. In the sum 
for the soft mode, po — 0, contributions from 2N — 1 
Matsubara modes are present for each loop. For the hard 
mode with po = 2ttT(N — 1), only N contributions are 
present, as \po + qo\ < 2tt(N—1)T. Thus different numbers 
of modes contribute and the counter-terms cannot be the 
same. This is always the case as long a.s N < oo. To surpass 
the problem in a constructive manner, each mode will be 
renormalized independently. 

Therefore, a counter-term Lagrangian is added, given 

by 

6C = 5m 2 A a (0) 2 + £ (5Z 3T (q )A^(q )A TfiU ( qo )A a v (q ) 



+ 8Z 3L (qo)A%(qo)A L ^(<lo)Al(qo) 
+ 6Z 3 (q )c a (q Q )d 2 c a (q )) 



1.25 



1.2 



1.15 — 



1.05 



1 — 

r i I 



n=0 mode 
n=1 mode 
n=2 mode 



_1_ 



10 



16 



18 



20 



Fig. 2. Dependence of Z% on the number of Matsubara fre- 
quencies N for the three lowest frequencies. Shown are the 
values for the k — 1/2 solution at T = 1.5 GeV. The solution 
for k ~ 0.4 is very similar, even on a quantitative level. 



the tree-level term is replaced 7 by l+SZ^^+Sm 2 jp 2 . Thus, 
it also generates a mass-renormalization. 

The last ingredient is the renormalization prescription, 
which will be taken as 



F(s) = 1, 



(16) 



where s is the subtraction point. Here s = T is chosen, see 
[25] for an alternative. This also ensures G(s) 2 Z(s) = 1 as 
is required in the 4-dimensional theory |21j . 

In case of H(0,p 2 ), two prescriptions are necessary. 
The first is 8 



lim \p\- 2 H(0.p 2 ) 



';;</ 



'9 



9 2 TC A r -^ 



(15) 



where A®(q ) are the modes of the gluon field and c a (qo) 
and c a (qo) are the modes of the ghost and anti-ghost field, 
respectively. A T j LpiV are the appropriate tensor structures 
of derivatives. 

This shortcoming is a consequence of the truncation of 
the Matsubara sum, yielding incorrect ultraviolet proper- 
ties: For large momenta the system under consideration is 
equivalent to a 3-dimensional theory of IN — 1 particles 
per 4d particle species. In such a theory all fields can be 
renormalized independently. When the limit N — > oo is 
performed, the wavefunction renormalizations will again 
coincide, as is shown 6 for Z 3 in Fig. [3 However, this pro- 
cess is logarithmically slow. 

In all equations, but the one for H(0,p), using the 
counter-terms of <|15ll amounts to replacing the tree-level 
term 1 by 1 + 6Z3 = Z3. In this way multiplicative renor- 
malization is obtained explicitly. In the equation for H (0, p) , 

6 For gluons, the value of the wavefunction renormalization, 
in the computationally accessible range of iV-values, is dom- 
inated by TV-dependent effects at mid-momenta. Hence, the 
logarithmic running is not easily discernible. 



(17) 



where 1713d is the tadpole-improved mass of the 3-dimensional 
theory 0. It depends on r = rrih/g 2 , which is again taken 
to be the same value as in ref. [S]. For fixed 53, this mass 
is independent of temperature. In addition, 



1 



H(0,s) 



1 



m 3d 

a 2 



(18) 



is required to fix the wavefunction renormalization. 

It should be noted that, by renormalizing at T, it 
is guaranteed that the correct 3-dimensional limit is ob- 
tained. This prescription requires that at T — > 00 all dress- 
ing functions approach 1 at infinity 0. 

Therefore, the explicit implementation of the renor- 
malization prescription for the DSE of a dressing function 
F with self-energy contributions / 



1 



Hp) 



1 + Hp) 



(19) 



7 Note that when using the projector JSJ, also Sm 2 is multi- 
plied by £. 

8 For numerical reasons, it is actually performed at 25 i, 
where Si is the numerical IR-cutoff of the integration. This 
causes only a marginal difference. 
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is then 



F(p) 
5Z 3 = -1(a) 

and in the case of H(0,p 2 ) 



= 1 + SZ 3 + I(p) 



H(0,p 2 ) 



1 + 5Z- 



I(P) 



Sm 2 = m 2 — lim p 2 I(p) 

p— >0 



SZ 3L = -1(a) + 



\\m. p ^, Q p 2 I(p) 



where m r = 7713d is the renormalized mass. 

This renormalization prescription amounts in total to 
a modified momentum subtraction scheme. 



6 Numerical results 
6.1 Propagators 

Here the full system of Eqs. (|7I9|I at all momenta will be 
treated 9 . The numerical method employed is discussed in 
ref. |2~51Ei2] . 

For N = 12 independent Matsubara frequencies, the 
result is shown in Fig. 01 for both sets of infrared expo- 
nents, k — 1/2 and K = 0.39760. There are several ob- 
servations. First of all, in the infrared the soft modes are 
nearly unaffected by the presence of the hard modes. The 
latter show a significant modification, compared to tree- 
level, although still only of the order of 30%. The plateau 
reached in the ultraviolet is an artifact of the truncation, 
yielding a finite wave-function renormalization. The hard- 
mode dressing functions exhibit some structure. There are 
maxima in all dressing functions at mid-momenta. This 
is most pronounced in the case of the dressing functions 
of the longitudinal gluon. These structures do not trans- 
late into a corresponding structure in the propagators, 
which are monotonically decreasing from a constant of or- 
der 1/pQ in the infrared to in the ultraviolet. For the 
hard-mode dressing functions it is also nearly irrelevant 
to which of the two soft infrared solution they belong. 

In general, the hard modes deviate less from unity 
with increasing n. The only exception seems to be the 
ultraviolet behavior of the n = 1 gluon modes. This is, 
however, an artifact of the truncation. If the Matsubara 
sum would not be truncated, all dressing functions would 
go to zero for sufficiently large momenta. On the other 
hand, the finite value attained due to the truncation can 



9 A separate discussion of a further truncation, the ghost- 
loop only truncation, can be found in ref. |25j . This truncation 
includes only tree-level and loop diagrams with at least one 
ghost-line. These results support that it is indeed the ghost 
sector which drives the infrared behavior of Yang-Mills the- 
ory, also at finite temperature, in accordance with the Gribov- 
Zwanziger scenario. 



be influenced by the requirement to fit continuously to 
the mid-momentum behavior. Hence, due to the peaks at 
mid- momenta in the gluon dressing functions, the ultra- 
violet plateau is increased for the n = 1 dressing func- 
tions, leading to the level reordering. If the peaks would 
be larger, this would also affect further modes. In the limit 
of N — > 00, the peaks reach a finite maximum and there- 
fore permit the plateaus all going to zero, thus restoring 
the correct behavior. 

Figure 01 displays the dependence of the full solutions 
on the number of Matsubara frequencies. While the soft- 
mode dressing functions are nearly unaffected, apart from 
the value of the renormalization constants, the effect on 
the hard mode dressing functions is significant. The ghost 
is quite insensitive, except for its wave-function renormal- 
ization. This is not the case for the gluons. The peaks 
at mid-momentum are sensitive to the number of Mat- 
subara frequencies included. The effect is largest for the 
longitudinal gluon dressing functions. From the available 
number of Matsubara frequencies, it is hard to estimate 
whether the peak grows to a finite value for an infinite 
number of Matsubara frequencies or not. It therefore can- 
not be excluded that the hard longitudinal gluon dressing 
functions violate the Gribov condition, once sufficiently 
many Matsubara frequencies are included. This would be 
very similar to the case of the transverse gluon dressing 
function in the infinite-temperature limit [5] and there- 
fore would necessitate a similar vertex construction for 
the longitudinal-gluon-transverse-gluon vertex as for the 
soft transverse gluon vertex. If such a vertex would be 
necessary, the result would be a finite peak, very similar 
to the present situation. Thus, the results would be even 
quantitatively quite similar. 

The disappearance of the peaks in the gluon dress- 
ing functions at N — 2 can be directly related to the 
vanishing of hard-mode couplings due to the restriction 
(po + 9o)/(27rT) < 2. Only hard-soft-mode couplings con- 
tribute. Thus, the peak is generated solely due to pure 
hard-hard interactions. 

In general, the n > 1 modes essentially follow the be- 
havior of the n = 1 mode, albeit much closer at tree-level. 

The dependence on temperature is shown in Fig. [5] 
The infrared coefficients for the soft modes of the ghost 
and transverse gluon dressing functions are affected, while 
the dressing function of the longitudinal gluon only changes 
slightly at mid-momenta. Therefore, the dressing function 
of the longitudinal gluon is dominated by its renormal- 
ized mass also when changing the temperature. As ex- 
pected, the hard-mode dressing functions become more 
and more tree-level like with increasing temperature. The 
peaks in the hard mode dressing functions shift to higher 
momenta, owing to the renormalization condition, and be- 
come smaller due to the increase of the effective mass po 
of the hard modes. At infinite temperature, they smoothly 
merge into the corresponding 3-dimcnsional solutions. At 
low temperatures, the opposite effect is observed. In gen- 
eral, the hard-mode dressing functions become more dy- 
namical as their effective mass decreases, albeit quite slowly. 
The insensitivity on the infrared solution of the soft sec- 
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Fig. 3. From top to bottom, the left panels show the dressing functions for the soft modes of the ghost and transverse gluon 
and the propagator of the longitudinal gluon. The right panels display the hard- mode dressing functions n — 1 to n — 11. 
Higher n correspond to lesser deviation from unity except for n = 1, see text. The solid lines are the ft = 1/2 solution and the 
dashed lines the n ^ 0.4 solution. Both are taken at T = 1.5 GeV, £ = £ = 1, and N = 12. Note the different momentum scale 
in the left and the right panels. 
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Fig. 4. The dependence of the solutions on JV at T = 1.5 GeV and £ = £ = 1. For the soft modes from top to bottom the left 
panels show the dressing functions of the ghost and transverse gluon and the propagator of the longitudinal gluon. The right 
panels give the dressing functions for the n — 1 hard mode. Only the « ~ 0.4 solution is shown. The n — 1/2 solution is similar 
and can be found in ref. |25j . The solid lines denote N — 2, dashed are N = 5, dotted N = 9 and dashed-dotted iV = 13. 
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Fig. 5. The dependence of the solutions on T at N — 5 and ^ = ^ = 1. The left panel shows the soft modes, where apart from 
the longitudinal gluon the dressing functions are shown. In the latter case the propagator is shown. The right panels show the 
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tor is not changed when reducing the temperature, and 
it may need a significantly lower temperature to induce a 
change. 



Employing t'Hooft scaling, the lowest temperatures 
that can be achieved with a numerically stable solution are 
of the order of 750 MeV. Using a temperature- independent 
renormalization scale and subtraction point |25) instead 
results merely in a shift in the onset of nonperturbative 
effects to larger momenta with some deformations, neces- 
sary to obey the renormalization conditions. In this case, 
it was also possible to cool down to T = 140 MeV, indi- 
cating the possibility of super-cooling. 

Varying the projection parameter £ leads to similar 
variations for the soft-mode dressing functions as in the 
3-dimensional case The largest effect is seen in the 
infrared. As the hard-mode dressing functions are insensi- 
tive to the infrared behavior of the soft modes due to their 
effective mass, there are only weak variations of them with 

Only at small £ an additional structure appears in the 
transverse gluon dressing function. This is likely due to the 
cross-term in Eq. ©. The effect of varying £ on the lon- 
gitudinal sector is negligible, as in the 3-dimensional case. 
Correspondingly, the effect of varying £ is only significant 
for the longitudinal sector. The sensitivity is much less 
pronounced then in the case of varying £. Especially, the 
soft longitudinal dressing function is nearly unaffected by 
variation of £. The reason is that the soft equation does 
not depend explicitly on £, as it can be divided out of 
the equation for £ ^ 0. Therefore, any £ dependence en- 
ters only indirectly by the weak dependence of the hard 
modes on £. Only in the case of £ ~ significant effects are 
found, owing to the pathology of the £ = case. Details 
of the dependence on £ and £ can be found in ref. |25j . 
The dependence on the three-gluon vertex construction 
is qualitatively not different from the infinite-temperature 
limit for the soft modes, and negligible for the hard modes 

EZ5I- 

At this point, a comparison to a different approach 
to obtain the high-temperature gauge propagators can 
be made. The usual continuum method is the semi-per- 
turbative hard thermal loop (HTL) approach It is 
based on resumming the hard mode contributions in self- 
energy diagrams. In the transverse infrared sector it is 
plagued by severe problems, due to its perturbative na- 
ture. The final result is a transverse gluon propagator 
with a particle like pole at p = 0. Thus the gluon ex- 
ponent t would vanish. This is in sharp contrast to the 
results found here and, as discussed in ref. such an in- 
frared behavior is not very likely. Also the lattice results 
in the infinite-temperature limit [5] and at temperatures 
somewhat above the critical 34 support the results found 
here. On the other hand, concerning the soft longitudinal 
mode and the hard modes, HTLs and the ansatz presented 
here find qualitatively similar results on the level of the 
propagators. 



6.2 Thermodynamic potential 

Concerning the contribution from the hard modes to the 
thermodynamic potential, it is problematic to use the ap- 
proach of Luttinger-Ward / Cornwall-Jackiw-Tomboulis 
(LW/CJT) [35], as done previously [51112] . In the present 
truncation scheme the 2PI contribution in the LW/CJT- 
action are neglected, and thus only the interaction part is 
considered. This latter part vanishes identically for a free 
system, i.e. for a system containing only tree-level dress- 
ing functions. This free contribution has thus to be added 
explicitly. Therefore, the hard-mode contribution, being 
essentially tree-level, can be calculated up to perturbative 
corrections to be the same as that of a non-interacting 
system of free gluons. The LW/CJT-expression does not 
capture the tree-level contribution of the soft modes, and 
it thus can be added here. This yields a thermodynamic 
potential of a gas of massless gluons with small correc- 
tions due to the explicit soft-mode contributions and the 
residual interactions of the hard modes. For N c = 3 it is 
given by pp 

G 1671-2 
lim — = — . (20) 



rp4 



90 



Using t'Hooft-scaling for the interaction strength, the cor- 
rections due to the soft contributions obtained by the 
LW/CJT action vanish most likely as 1/T |5H25| . The only 
remaining contribution in the infinite-temperature limit is 
hence the Stefan-Boltzmann-contribution l|20|l . in agree- 
ment with results from lattice calculations p^. Therefore, 
the thermodynamic properties of the high-temperature 
limit are governed by the hard modes. However, as the 
calculation of the thermodynamic potential is problem- 
atic in the current approach [R*ll25| . this can only be taken 
as an indication for such a behavior. 

At temperatures of the order of the phase transition, 
the contribution of the soft modes to the thermodynamic 
potential is probably non-negligible. Furthermore, as dis- 
cussed in ref. (21], these contributions can be highly rele- 
vant to the pressure p and lead to a non-vanishing value 
for the trace of the stress-energy tensor, e — 3p, where e 
is the energy density. The vanishing of the trace would be 
expected for an ideal gas. Such an effect cannot be pro- 
vided by perturbation theory |36|. but has been observed 
in lattice calculations 



6.3 Schwinger functions 

As stated in section^ an unambiguous sufficient (although 
not necessary) signal for the confinement of a particle is a 
violation of positivity. This can be investigated using the 
Schwinger function, defined as 



A(z) = 



1 



dp cos( zp 3 



2 " 

n 



(21) 



where F is either Z or H. In the infinite-temperature limit 
a clear signal for positivity violation has been found for 
the transverse gluon [SJ, in agreement with lattice results 
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Fig. 6. The Schwinger function for the gluon propagator. The left panels show the results for the soft modes and the right 
panels for the hard modes with n = 1 and n = 2. The top panels give the results for the transverse propagator and the bottom 
panels for the longitudinal propagator. The temperature is T = 1.5 GeV, f = £ = 1, and N — 5. The result has been obtained 
with roughly half a million Fourier frequencies 1251 . 



|51llUII37| . However, for the longitudinal gluon no unam- 
biguous signal could be found. It is therefore interesting 
to repeat the study at lower temperatures. The results are 
shown in Fig. [B] 

The soft transverse function does show a clear sig- 
nal of positivity violation. The hard modes do not show 
a clear sign of oscillation within the available precision. 
Compared to the high-temperature much clearer 

signal for positivity violation is found in the soft longitudi- 
nal propagator. Furthermore, the longitudinal Schwinger 



function is much more similar to the transverse one. There- 
fore, the contribution from the hard modes seem to be 
negative, shadowing the oscillations induced by the soft 
modes in the infinite-temperature limit. It is also this 
observation, which restricts the reliability of the result. 
Nonetheless, these results substantiate the non-triviality 
also of the longitudinal sector and requires further atten- 
tion. This is especially true for the soft longitudinal mode, 
where it is yet unclear whether the positivity violation is 
an artifact of the truncation and if not what is its origin. 



Axel Maas et al.: The High- Temperature Phase of Landau-Gauge Yang- Mills Theory 



13 



Known is so far that the extent of such violations depend 
on the interaction strength [51l25| . In the case of the trans- 
verse gluon, the violation is driven by the vanishing of the 
propagator in the infrared, and thus a quite robust state- 
ment. 



7 Discussion 

Before discussing further the implications of the results 
found here on the structure of the Yang-Mills phase dia- 
gram, it is worthwhile to take a step back to gather and 
assess the findings. Especially, the reliability of the results 
deserves special attention. 

The reliability of the exploratory study of finite tem- 
perature effects presented here and thus of the proper- 
ties of the hard modes is a major concern. The success 
of the truncation scheme in the infinite-temperature limit 
is based on its exactness in the ultraviolet and its pre- 
sumed exactness in the infrared. Thus the consequences 
of the deficiencies at intermediate momenta are strongly 
constrained. This constraint is lost by cutting off the Mat- 
subara sums. As a consequence, the ultraviolet properties 
of the hard modes are incorrect and a rooting in perturba- 
tion theory is prevented. This amplifies the second prob- 
lem. The hard modes do not reach into the 4-dimensional 
infrared region due to their effective mass. Therefore the 
advantages of the truncation scheme are lost to a large 
extent and the equations for the hard modes become un- 
reliable. 

Nonetheless, the presented results are relevant. At suf- 
ficiently large temperatures, the truncation-dependent self- 
energies of the hard modes are suppressed by the large 
effective mass and they are dominated by the truncation- 
independent tree-level terms. Thus, the truncation arti- 
facts vanish when the temperature goes to infinity. The 
difficult part is then to assess down to which temperature 
the results are still reliable, at least qualitatively. 

As it is found that the system is only very weakly de- 
pendent on temperature, qualitative conclusions can prob- 
ably be drawn even at quite small temperatures. This is 
due to the fact that the infrared properties of the soft 
modes are nearly independent of the hard modes. In gen- 
eral, the hard modes effectively decouple because of their 
effective mass, which is large compared to Aqcd, even at 
the phase transition. The only exception is the genera- 
tion of the screening mass for the soft longitudinal gluon, 
which can be found on quite general grounds and thus can 
be considered reliable. This is also supported by the sys- 
tematic error estimations performed in refs. pll25| . which 
do not find qualitative effects but only small quantitative 
ones. 

Although the quantitative results may be subject to 
change, the qualitative result is expected to be quite re- 
liable. Probably down to temperatures within the same 
order of magnitude as the phase transition, the high tem- 
perature phase consists of strongly interacting soft modes 
exhibiting confining properties and inert hard modes. This 
is supported by lattice results, which find that the infinite- 
temperature limit is effectively reached at about 2T C [H] 



and the qualitative features of the infinite-temperature 
limit up to the highest temperatures, for which propaga- 
tors are available, of 6T C This constitutes the major 
result found in this work. 

Concerning the thermodynamic potential, no final con- 
clusion can be drawn. The results indicate that the gross 
thermodynamic properties far away from the phase transi- 
tion are completely dominated by the hard modes and that 
a Stcfan-Boltzmann behavior is reached in the infinite- 
temperature limit. Considering the crude assumptions 
made and the difficulties encountered, this result is in- 
dicative at best. Nonetheless, if it can be substantiated, 
this would allow to understand how a Stefan-Boltzmann- 
like behavior can emerge from a non-trivial theory. At 
the same time, the results also indicate that, at low tem- 
peratures, nonperturbative effects will play a role in the 
thermodynamics and certainly will especially be relevant 
in the vicinity of the phase transition. 

8 Concluding remarks 

Looking at the results found in ref. [51l25| and the solu- 
tions at zero [2] and small temperatures ^21> a coherent 
picture emerges. Still, this picture is overshadowed by the 
problems of truncation artifacts, especially at finite tem- 
peratures. 

The main difference between the low-temperature and 
the high-temperature phase is not primarily one between 
a strongly interacting and confining system and one with 
only quasi- free particles. The chromoelectric gluons, whose 
infrared behavior change from over-screening to screen- 
ing, come somewhat close to such a picture, and the hard 
modes certainly do. The latter become free due to the vac- 
uum property of asymptotic freedom and the dependence 
of their energy on the temperature, thus in an expected al- 
though not entirely trivial manner. The chromomagnetic 
gluons remain over-screened in the infrared and are thus 
confined. In the sector transverse to the heat-bath, the 
results provide strong evidence for the Zwanziger-Gribov 
and/or the Kugo-Ojima scenarios, even at very high tem- 
peratures. 

Including the observation of super-heating in the low- 
temperature regime |12) and super-cooling at high tem- 
peratures, the scenario emerging is a chromoelectric phase 
transition of first order. This phase transition connects two 
different strongly interacting phases of Yang- Mills theory. 
Especially quantum fluctuations always dominate thermal 
fluctuations and at least part of the gluon spectrum is al- 
ways confined. 

In the vicinity of the phase transition, the consequences 
of nonperturbative effects are likely also relevant to ther- 
modynamic properties, especially to the pressure and the 
trace anomaly |36j . This underlines the importance of non- 
perturbative effects at least for the temperature range rel- 
evant to experiment. 

In summary, evidence is found that, in the high-temper- 
ature phase of Yang-Mills theory, the elementary excita- 
tions exhibit quite non-trivial correlations. The soft in- 
teractions are non-trivial at all temperatures, and the de- 
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coupling of the hard modes is evident. This establishes the 
main features of the high-temperature phase of Yang-Mills 
theories in Landau gauge. The results found here com- 
ply with Gribov-Zwanziger or Kugo-Ojima type confine- 
ment mechanisms for the dimensionally reduced theory. 
Therefore, the high-temperature phase is definitely non- 
perturbative. It will be interesting to investigate the con- 
sequences for quarks and their correlations at high tem- 
perature on the one hand and for cosmology and the early 
universe on the other. 

The authors are grateful to Burghard Griiter for valuable dis- 
cussions and to Jan M. Pawlowski for a critical reading of the 
manuscript and helpful remarks. This work was supported by 
the BMBF under grant numbers 06DA116 and by the Helmholtz 
association (Virtual Theory Institute VH-VI-041). 



A Running coupling 

For fixing the coupling there are two options to consider. 

If fj. is fixed, gs = g 2 T grows without limit as T grows. This 
does not necessarily pose a problem, as the infinite-temperature 
propagators are independent of gs as long as they are expressed 
as a function of p/g 2 0. Under such circumstances, all mo- 
menta below T are effectively infrared, and the nonperturba- 
tive regime would extend to all momenta. Here this possibility 
will not be pursued any further 10 . 

Alternatively, it is possible to use the limiting prescrip- 
tion g 2 (fi(T))T = Coo, where Coo is an arbitrary constant, ef- 
fectively performing a renormalization group transformation 
when changing the temperature. This fixes g 2 to be propor- 
tional to Aqcd, the dynamical scale of QCD |31|. while the 4- 
dimensional coupling vanishes like 1/T. This defines a t'Hooft- 
likc scaling in T. Note that g — *• for T — > oo corresponds to 
the conventional arguments of a vanishing coupling in the high- 
temperature phase |33) . However, as in the large-iV c limit, a 
t'Hooft-scaling is performed, generating a well-defined theory. 
Hence, this possibility defines a smooth 3-dimensional limit 
with a finite 3-dimensional coupling constant. Here Coo = 1 is 
chosen. 

Furthermore, it would be useful to give explicit units for 
the temperature scale. In the vacuum, the scale is fixed via 
a comparison of the running coupling to perturbation theory 
|23) . This is not possible here because the Matsubara sum is 
truncated, and the coupling cannot be calculated reliably for 
momenta of the order of 2-nNT. 

The simplest procedure is to compare gz to lattice calcu- 
lations. There, the effective 3-dimensional coupling is found to 
be g 2 (2T c ) /2T C = 2.83 0. The phase transition temperature is 
T c = 269±1 MeV [T|. Albeit this will not be exactly g 2 -2T c due 
to the truncation, a first approximation is to require the same 
temperature scale at 2T C . Using the fixed g 2 T = Coo prescrip- 
tion, then 1 = Coo/(2T c • 2.83). Since here Coo = 1 in internal 
units, internal units have to be multiplied by 1.5 GeV to yield 
physical units. This temperature scale is used in section^] The 
ratios of temperatures are independent of this prescription. 



See |25j for a thorough discussion of this case 



B Subtraction of spurious divergences 

The spurious divergences in contributions of hard modes due 
to the integration are dealt with in exactly the same way as in 3 
dimensions by adjusting the tadpole terms. Each integration 
kernel K in the longitudinal equation is split as 

K = Ko + K D (22) 

where Ko is finite and Kd divergent upon integration. Each 
kernel Kd will be compensated by the corresponding tadpole 
in 

A new kind of spurious divergences appear when perform- 
ing the Matsubara sum. When including more and more Mat- 
subara frequencies, it is found that their contribution in the 
gluon equations ® and Q at \p\ < max(qo) scales as qo, thus 
behaving as a quadratic divergence. This is an artifact of using 
a finite number of Matsubara frequencies and does not van- 
ish as long as the number of frequencies is finite, no matter 
how large the number. This behavior is spurious and must be 
removed. 

This can be performed by the replacement 

F{q, qo)F(p + q, q + pa) 
-> (F(q,qo) - (F(q + p,qo+p ) - , (23) 

in all affected loops in the gluon equations. F is a generic dress- 
ing function and Z% its wave-function renormalization. 

In the transverse equations all spurious divergences, includ- 
ing the one of the Matsubara sum, are again removed at £ = 3. 
Therefore, the replacement 112MB will only be necessary in those 
contributions which are proportional to (£ — 3). Hence the in- 
tegral kernels are split differently than in 122H as 

K(0 = K + (t-3)K 3 + Ko(0- (24) 

Here Ko and Kz are finite and independent of £, and Kd con- 
tains all divergences upon integration. The corresponding sub- 
traction is then performed by the replacement 

KF(q)F(p + q)^K F(q)F(p + q) (25) 

Ht-3)K 3 (F { q)-±-) (F(p + q) -±). 

In the soft transverse equation it is necessary in addition 
to remove tadpole- like structures at £ 7^ 3 in K3 |25) . which 
alters 125H to 

KF(q)F(p + q)^ K F(q)F(p + q) 

+(C-3) (^-^(lim^oP 2 ^)) • 
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